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The recent WMAP data have confirmed that exotic dark matter together with the vacuum energy 
(cosmological constant) dominate in the flat Universe. The nature of the dark matter constituents 
cannot be determined till they are directly detected. Recent developments in particle physics provide 
0^ ' a number of candidates as constituents of dark matter, called Weakly Interacting Massive Particles 

\ (WIMPs). Since these interact weakly and are of low energy, they cannot excite the target and can 

only be detected via measuring the recoiling nucleus. For all WIMPs, including the most popular 
candidate, the lightest supersymmetric particle (LSP), the relevant cross sections arise out of the 
' following mechanisms: i) The coherent mode, due to the scalar interaction, ii) The charge coherent 

D , mode, with only proton contribution, as in the recent case of secluded dark matter scenario and 

iii) The spin contribution arising from the axial current. In this paper we will focus on the spin 
contribution, which maybe important, especially for light targets. 
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^ . cold dark matter (CDM): 



INTRODUCTION 



Combining the recent WMAP data [l[ with other experiments one finds that most of the matter in the Universe is 



VL b = 0.0456 ± 0.0015, Qcdm = 0.228 ± 0.013, Ct A = 0.726 ± 0.015 



Since the non exotic component cannot exceed 40% of the CDM there is room for exotic WIMPs (Weakly 
Interacting Massive Particles). In fact the DAMA experiment Q has claimed the observation of one signal in direct 
detection of a WIMP, which with better statistics has subsequently been interpreted as a modulation signal [4| . These 
data, however,if they are due to the coherent process, are not consistent with other recent experiments, see e.g. 
EDELWEISS (I, CDMS @ and XENON 0). The DAMA results could still be interpreted as due to the spin cross 
, section, but with a new interpretation of the extracted nucleon cross section. 

Supersymmetry naturally provides candidates for the dark matter constituents (see, e.g, the review jgj). In the 
most favored scenario of supersymmetry the LSP (Lightest Supersymmetric Particle) can be simply described as a 
Majorana fermion, a linear combination of the neutral components of the gauginos and higgsinos @],@. In most 
calculations the neutralino is assumed to be primarily a gaugino, usually a bino. Models which predict a substantial 
fraction of higgsino lead to a relatively large spin induced cross section due to the Z-exchange. Such models have been 
less popular, since they tend to violate the relic abundance constraint. These fairly stringent constrains, however, 
apply only in the thermal production mechanism. Furthermore they do not affect the WIMP density in our vicinity 
derived from the rotational curves. Thus one may assume that large spin cross sections are possible in such models 
0, [ll|, which are non-universal gaugino mass models and give rise to large higgsino components. Sizable spin cross 
sections also arise in the context of other models, which have appeared recently [I3|-jl3| (see also Ellis et al 15 1 



for an update). 

Spin induced cross sections also arize in the case Kaluza-Klein (K-K) WIMPs in models with Universal Extra Dimen- 
sions (UED) [lfjj]. This occurs regardless of whether the WIMP is a K-K boson or a K-K neutrino. They can also arise 
in technicolor theories [l7| • In the Ultra Minimal Walking Technicolor model ID, 19|, 2(| there exist singlet composite 



Majorana fcrmionic states. These, taken as dark matter candidates, can lead to spin induced cross sections. 

Knowledge of the spin induced nucleon cross section is very important since, for some special targets, it may lead 
to transitions to excited nuclear states, which provide the attractive signature of detecting the de-excitation 7 rays 



in or without coincidence with the recoiling nucleus |21|. 22j. Furthermore it may dominate in light systems like 3 He 
and 19 F, which offer some attractive advantages 23 1. [241]. 



In light of the above it is clear that the spin mechanism needs be considered. In this article we will discuss the 
theoretical ingredients needed to obtain the WIMP-nuclear spin induced cross sections. Then we will give expressions 



for and calculate the event rates, both modulated and unmodulated, in terms of the elementary proton (a p ) or neutron 
(<r n ) cross sections. After that we will provide exlusion plots in the (<t p , <7 n ) plane in terms of parameters relevant to 
the experiments, for various targets of experimental interest. Our results will be presented in a way that will make 
them useful in the analysis of the data of the odd mass targets. 



THE ESSENTIAL THEORETICAL INGREDIENTS OF DIRECT DETECTION. 



Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational 
curves, it is essential to directly detect such matter. Such a direct detection, among other things, may also unravel 
the nature of the constituent of cold dark matter, namely the Weakly Interacting Massive Particle (WIMP). The 
possibility of such detection, however, depends on the nature of its constituents. Our main conclusions apply to all 
heavy WIMPs. Since the WIMP is expected to be very massive, m x > 30GeV, and extremely non relativistic with 
average kinetic energy T as 50KeV(m x /l00GeV), it can be directly detected mainly via the recoiling of a nucleus 
(A,Z) in elastic scattering. The event rate for such a process can be computed from the following ingredients: 

1 . An effective Lagrangian at the elementary particle (quark) level obtained in the framework of supersymmetry 
as described , e.g., in Refs @, H^j- An analogous procedure can be found in the case of K-K WIMPs in Universal 
Extra Dimension (UED) models and technicolor theories [13] • 

2. A well defined procedure for transforming the amplitude obtained using the above mentioned effective Lagrangian 
from the quark to the nucleon level where This step is not trivial, since the obtained results depend crucially 
on the content of the nucleon in quarks other than u and d. This is particularly true for the scalar couplings, 
which are proportional to the quark masses Hi, [53, as well as the isoscalar axial coupling. 



3. 



4. 



Nuclear matrix elements. 

These must be obtained with as reliable as possible many body nuclear wave functions. Fortunately in the most 
studied case of the scalar coupling the situation is quite simple, since then one needs only the nuclear form 
factor. Some progress has also been made in obtaining reliable static spin matrix elements and spin response 
functions 

A velocity distribution for WIMPs 

In this article we will follow the standard practice and assume a M-B distribution, but other perhaps more 
realistic velocity distributions have also recently been considered 3^, 34 1 



Since the obtained rates are very low, one would like to be able to exploit the modulation of the event rates due to 
the earth's revolution around the sun M 0, Ezl, M 



THE WIMP NUCLEUS CROSS SECTIONS 



The standard (non directional) differential rate can be written as 

dR=^-^-da(u,v)\v\, (1) 

where m is the detector mass, p(0) = 0.3GeV/cm 3 is the WIMP density in our vicinity, v its velocity and m x o its 
mass and da(u,v) is given by 

da(u,v) — 2 ^ U bv)2 [ZsF 2 (u) + t spm F n (u)] , (2) 

where u is the energy transfer Q in dimensionless units given by 

U =TT > Qo= [m p Ab}- 2 = 40/U 4/3 MeV, (3) 

with b is the nuclear (harmonic oscillator) size parameter. F(u) is the nuclear form factor and Fn(u) is the spin 
response function associated with the isovector channel. The scalar cross section is given by: 
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FIG. 1: The spin response functions Fu(Q), F i(Q) and F 00 (Q) in the case of the target 9 F as a function of the 
energy transfer. In the region of interest for dark matter searches they are not distinguishable. 



&fj v o is the WIMP-nucleon scalar cross section. Note that, since the heavy quarks dominate, the isovector contribution 



is negligible, i.e. the proton and nucleoli cross sections are the same. The spin Cross section is given by: 
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The spin response functions Fij , properly normalized to unity at momentum transfer zero, in the energy transfers of 
interest are almost the same for all isospin channels = 0, 1 . As a matter of fact in the case of 19 F we get: 



Fu{u) 



0.0119u 4 - 0.1450u J + 0.6620?r - 



Au 



1 



(7) 



Fqi = e 



0.0124u 4 -0.1487iT 



o 4it 
0.6677u 2 - — 



(8) 



F 00 = e~ u ( 0.0131u 4 - 0.1525u 3 + 0.6733u 2 - — + 1 



(9) 



Thus they are indistinguishable (see Fig. [T]) for the energy transfers of interest. Hence 



S(u) « s(o) = ^Oo(o) + ni(o)) 



(10) 



The nuclear matrix elements r2i (0) (f2o(0)) associated with the isovector (isoscalar) components are normalized so 



that, in the case of the proton, they yield 



spin 



1 at u = 0. 



The couplings f A {f A ) are obtained by multiplying the corresponding elementary amplitudes obtained at the quark 
level by suitable renormalization factors g A and g\ given in terms of the quantities Aq given by Ellis [39j 



g° A = Au + Ad + As = 0.77 - 0.49 - 0.15 = 0.13 , g\ = Au - Ad = 1.26 



(11) 



Thus, barring very unusual circumstances at the quark level, the isovector component is expected to be dominant. It 
is for this reason that we started our discussion in the isospin basis. 

Heavy nuclei, however, are theoretically described in terms of protons and neutrons and the experiments are also 
analyzed this way. So will present our results in this basis. The proton and neutron cross section are given by: 



<$ = 3a \f° A + f\\ 2 = 3<roK\ 2 , <5 = 3a \f A - f\\ 2 = 3a \a n \'< 



(12) 
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TABLE I: The static spin matrix elements for various nuclei. For 3 He see Moulin, Mayet and Santos [41]. For the 
other light nuclei the calculations are from DIVARI j^]. For 73 Ge and 127 I the results presented are from Ressel et 
al [111 (*) and the Finish group et al [HI (**). For 207 Pb they were obtained by the Ioannina team (+). [H, [13]. 





3 He 


19 p 


29 Si 


23 Na 


73 Ge 


127j* 


127j** 


207p b + 


fio(O) 


1.244 


1.616 


0.455 


0.691 


1.075 


1.815 


1.220 


0.552 


fii(O) 


-1.527 


1.675 


-0.461 


0.588 


-1.003 


1.105 


1.230 


-0.480 


fip(O) 


-0.141 


1.646 


-0.003 0.640 


0.036 


1.460 


1.225 


0.036 


n„(o) 


1.386 


-0.030 


0.459 


0.051 


1.040 


0.355 


-0.005 


0.516 






2.91 


-0.50 


2.22 














2.62 


-0.56 


2.22 










Vth( s P in ) 




0.91 


0.99 


0.57 











with a p and a n are the proton and neutron spin amplitudes, which, of course, depend on the model. In the case of 
the LSP [H 



(Jq = -!-(G F m„) 2 = 0.77 x l(T 38 cm 2 = 0.77 x 10" 2 pb. 
In extracting limits on the nucleon cross sections from the data we will find it convenient to write: 

^ r \2_svin _svin 1 i<~i ^— , o ^— — „ i <5 1 2 1 i 



s spm = (-^) 2 KT > <r = -\n p ^; + n n ^ 5 \ 2 = -\\n p \^+ \n n \^ e ^+ s ^\ 2 , (13) 

f-r(P) 3 3 

where O>(0) and 0(0) are the proton and neutron components of the static spin nuclear matrix elements, 8a is the 
relative phase between them (zero or tt) and 5 the relative phase between the amplitudes a p and a n . 
The nuclear spin ME are defined as follows: 



0,d») = -< J J\<?z(p)\J J >- > 0(0) = J —J— ~< J J\&z{n)\J J >- (14) 



where J is the total angular momentum of the nucleus and <j z = 2S Z . The spin operator is defined by S z (p) = 
^ i=1 iS z (i), i.e. a sum over all protons in the nucleus, and S z (n) — Yli=i i.e. a sum over all neutrons. 

Furthermore 

fio(0) = fi p (0) + 0,(0) , Oi(0) = O p (0)-0(0). (15) 

The spin ME can be obtained in the context of a given nuclear model. Some such matrix elements of interest to the 
planned experiments are given in table |H The shown results are obtained from DIVARI j2{|, Ressel et al (*) [13], the 
Finish group (**) [3l[ and the Ioannina team (+) |4(J, [32| . 

Before concluding this section we should emphasize that from the spin matrix elements of Table Q] those associated 
with 19 F are the most reliable for the following reasons j29|: 

• The light s-d nuclei are very well described within the interacting shell model. 

• The magnetic moment of the ground state is dominated by the spin (the orbital part is negligible). 

• The calculated magnetic moment is quite large and in good agreement with experiment. 

To summarize: The proton and neutron spin cross sections can be obtained in a given particle model for the WIMP's. 
As we have seen there is a plethora of such models to motivate the experiments. Some of them may yield as high as 
a few tens of events per kg of target per year [12] ■ But most of them depend on imput parameters that are not well 
detemined. So none of them seems to be universally accepted. Thus in the present work, rather than following the 
standard procedure of providing constrained parameter spaces, we will treat the proton and neutron cross sections as 
parameters to be extracted from the data. This can be done, once the nuclear spin matrix elements are known, for 
various values of the phase difference S. The only particle parameter we will retain is the WIMP mass, which is the 
most important, since it enters not only in the elementary cross sections but the kinematics as well. 
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FIG. 2: On the left the quantity c sp i n (A, ji r (A), m x o) £ S pm for the A=19 system is shown for two cut off values 
Qmin = 0, continuous curve, and Q m in = 10 keV, dotted curve. On the right the same quantity is shown for the 
A=127 system. The advantages of the lighter target, especially for light WIMP, are obvious. 



EXPRESSIONS FOR THE RATES AND SOME RESULTS 



To obtain the total rates one must fold the diffrential rate of Eq. {T]) with WIMP velocity and then integrate the 
resultin expression over the energy transfer from Q m in determined by the detector energy cutoff to Q ma x determined 
by the maximum WIMP velocity (escape velocity, put in by hand in the M-B distribution), i.e. v esc = 2.84 vq with 
vq the velocity of the sun around the center of the galaxy (229 Km/ s). 

Ignoring the motion of the Earth the total (non directional) rate is given by 



R = Rt(a, Q mm ) , R=PW™ ( J* ^fl) A 2 + ^ 



(16) 



The WIMP parameters have been absorbed in R. The parameter t takes care of the nuclear form factor and the 
folding with WIMP velocity distribution l^, H, 52] (for its values see table ITT|) . It depends on Q mm , i.e. the energy 
transfer cutoff imposed by the detector and a = \iirbvoV2\~ 1 ■ 
In the present work we find it convenient to re-write it as: 



R = K(a 1 ) 



^COfl 

(A, /j, r (A), m x o 



, a P,x° 
01 



t (A,/i r (yl),m x o) 



-spin 



(71 



For the spin cross section it is convenient to take a\ = 10 5 pb. Thus 

P(0) 



K{a 1 ) 



100 GeV : 



1.60 10" 2 y- 1 



P(0) 



V¥ 2 



0.3GeVcm~ 3 IKg 280fcms" 



(17) 



(18) 



For the coherent mode it may be more convenient to pick a\ = 10 7 pb, which is close to the present experimental 
limit. Furthermore 



c coh (A,fi r (A),m, x a)) 



100 GeV 



Hr(A) 



Hr{p) 



A t coh {A) , C S pin (A, fj, r (j4) , m x o)) = 



lOOGeV 



t 



sp/n 



(A) 



A 



(19) 



The parameters c co ) l {A 1 fi r (A), to x o), c sp i n {A 1 fi r (A), to x o), which give the relative merit for the coherent and the spin 
contributions in the case of a nuclear target compared to those of the proton, are tabulated in table |TT] for energy 
cutoff Q m in = 0, 10 keV. Thus via Eq. (fl~7|) we can extract the nucleon cross section from the data. 

The quantity c spm (A, /Lt r (A), m x o) Cspin, when the isoscalar contribution is neglected and employing O 2 = 1.22 (2.81) 
for 127 1 ( 19 F), is shown in FigH In the case of the spin induced cross section, the light nucleus 19 F has certainly an 
advantage over the heavier nucleus 127 I (see Fig. [2]). For the coherent process, however, the light nucleus is disfavored, 
(see Table HJ. 
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TABLE II: The factors cl9 = c co /j(19, /x r (19), m x o), sl9 = c sp i n (19, /i r (19), m x o), cl9 = c co /i(73, /x r (73), m x ), 
s73 = c sp i„(73, /Lt r (73), m x ) and cl27 = c co ? l (127, /it r (127), m x o), sl27 = c sp i„(127, jU r (127), m x o) for two values of 
Qmin- Also given are the factors s3 = c S pi n (3, £t r (3), m x ) for Q mm = 0. 



Qmin 




m x (GeV) 


keV 




20 


30 


40 


50 


60 


80 


100 


200 





t(3,s) 


1.166 


1.166 


1.166 


1.166 


1.166 


1.166 


1.166 


1.166 





c3 


131 


92.6 


71.6 


58.3 


49.2 


37.5 


30.3 


15.4 





s3 


14.6 


10.3 


7.95 


6.48 


5.47 


4.16 


3.36 


1.71 





t(19,c) 


1.153 


1.145 


1.138 


1.134 


1.130 


1.124 


1.121 


1.112 





t(19,s) 


1.132 


1.117 


1.105 


1.096 


1.089 


1.079 


1.072 


1.056 





cl9 


11500 


10500 


9420 


8500 


7700 


6450 


5540 


3210 





sl9 


31.2 


28.3 


25.4 


22.8 


20.6 


17.2 


14.6 


8.40 





t(23,c) 


1.107 


1.099 


1.092 


1.089 


1.085 


1.079 


1.076 


1.068 





t(23,s) 


1.075 


1.061 


1.050 


1.041 


1.035 


1.025 


1.018 


1.003 





c23 


16100 


15200 


14100 


1300 


11900 


10200 


8830 


5280 





s23 


29.5 


27.8 


25.6 


23.4 


21.4 


18.2 


13.8 


9.45 





t(73,c) 


1.119 


1.083 


1.047 


1.014 


0.984 


0.933 


0.893 


0.780 





t(73,s) 


1.135 


1.112 


1.088 


1.064 


1.043 


1.006 


0.976 


0.886 





c73 


113000 


131000 


139000 


143000 


142000 


187000 


130000 


935000 





s73 


20.8 


23.9 


25.2 


25.5 


25.2 


23.9 


22.2 


15.4 





t(127,c) 


0.984 


0.844 


0.721 


0.621 


0.542 


0.430 


0.358 


0.213 





t(127,s) 


0.948 


0.796 


0.671 


0.574 


0.501 


0.401 


0.340 


0.220 





cl27 


206000 


225000 


223000 


211000 


197000 


169000 


145000 


82400 





sl27 


12.3 


13.1 


12.8 


12.1 


11.3 


9.7 


8.5 


5.3 


10 


t(19,c) 


0.352 


0.511 


0.592 


0.639 


0.667 


0.710 


0.720 


0.773 


10 


t(19,s) 


0.340 


0.489 


0.563 


0.606 


0.631 


0.669 


0.676 


0.720 


10 


cl9 


3500 


4676 


4902 


4789 


4546 


4075 


3557 


2233 


10 


sl9 


9.3 


12.4 


12.9 


12.6 


11.9 


10.6 


9.3 


5.8 


10 


t(73,c) 





0.020 


0.119 


0.246 


0.363 


0.539 


0.651 


0.847 


10 


t(73,s) 





0.0175 


0.105 


0.213 


0.311 


0.453 


0.539 


0.677 


10 


c73) 





2310 


15300 


32900 


49600 


73400 


86300 


89300 


10 


s73 





0.39 


2.5 


5.3 


7.9 


11.6 


13.4 


13.4 


10 


t(127,c) 


0.000 


0.156 


0.205 


0.222 


0.216 


0.191 


0.175 


0.109 


10 


t(127,s) 


0.000 


0.135 


0.177 


0.192 


0.190 


0.174 


0.165 


0.121 


10 


cl27 





41500 


63200 


75500 


78500 


74900 


71000 


42200 


10 


sl27 


0. 


2.2 


3.4 


4.0 


4.3 


4.2 


4.1 


2.9 



TABLE III: The experimental sensitivity ratios (ESR) for various targets assuming a WIMP mass of 50 GeV. p, n 
and iv correspond to the elementary proton, neutron and isovector dominance respectively. 





3 He 


19 p 


23 Na 


73 Ge 


127j 


p 


1.4xl0 3 


4.1xl0 2 


2.7xl0 3 


1.3xl0 7 


2.5xl0 4 


n 


1.4x10 


1.2xl0 b 


4.3x10° 


1.6xl0 4 


4.2x10° 


IV 


1.1x10 


4.0xl0 2 


3.2xl0 3 


1.7xl0 4 


4.3 xlO 4 



The experimental sensitivity ratios (ESR), i.e. the extracted from experiment nucleon cross section ratios satisfy: 



ESR = 



'k,x° 



pa" 



oh(A,Vr(A),m x o) 

tin 

(A, /j, r (A),m x a) 



A k 

n 9 * ' 



p, n, iv, for proton, neutrom, isovector respectively (20) 



The quantity ESR for a WIMP mass of 50 GeV is shown in table IIII1 It is clear from this table why the limits on 
the spin cross section extracted from all targets is much bigger compared to that extracted for the coherent mode. 



G 



We should emphasize that the elementary cross sections do not depend on the target. It is only the values extracted 
from experiment that do so, giving a measure of the sensitivity of the various experiments. The elementary cross 
sections only depend on the particle model and the structure of the nucleon. Thus, e.g., in the case of K-K WIMPs 
the coherent cross section dominates, if the WIMP is a K-K gauge boson, but the spin cross section is bigger, when 
the WIMP is a K-K neutrino 

If the effects of the motion of the Earth around the sun are included, the total non directional rate is given by 



R = K{ax) 



c coh (A, fJLr(A), m x o) — ; — (1 + h(a, Q mm )cosa) 



(coherent), 



(21) 



R = K(ax) 



l - Spl II 

(A, n r (A),m x o) 



nuc 



(1 



i (a, 



)cosa) 



(spin), 



(22) 



where h (h sp i n ) are the modulation amplitudes and a is the phase of the Earth, which is zero around June 2nd. We 
are going to only briefly discuss the modulation amplitudes here since they depend only on the WIMP mass and 
are independent of the other particle parameters. In the case of the two very light targets, however, they are pretty 
independent of the WIMP mass. In fact for the light systems : 



h = h spm = 0.0232 (for A=3) and h = 0.0229, h s 



0.0227 (for A=19). 



(23) 



Actually for the A=19 system there is about (10%) reduction as the WIMP mass increases. 

In the case of the target 3 He the quantity t is also essentially independent of the WIMP mass, since the WIMP is 
expected to be much heavier than the nuclear mass. From table |TT] we see that the coherent rate is quite small for 
this light system, but the spin induced rate is only a factor of two smaller than that for 19 F. As we have already 
mentioned this nucleus, has definite experimental advantages fill ]. 

In many instances the experiments are interested in the differential event rate. This is a function of two variables, 
the WIMP mass and the energy transfer Q. For the light systems, however, the dependence on the WIMP mass is 
rather weak, especially for heavy WIMPs. Thus the presentation of the results is relatively simple and we are going 
to present them here. One finds: 



dR irt \ 



dR Q (Q,A,m^ 
dQ 



(1 + H (Q, A, m x ) cosa) 



01 



(24) 



with an analogous expression for the coherent mode. The time average quantity dRo{Q , A, m x ) / dQ and the relative 
modulation amplitude H(Q, A,m x ) are shown in Fig. [3j The differential cross section is normalized so that the area 



under the corresponding curve gives the value c sp i n (A, fi r (A), m x o ) of Eq. (|22p . Note that the quantity H, being the 
ratio of two amplitudes, the amplitude for modulation divided by the time independent amplitude, is independent of 
the nuclear model. So it is the same for the spin and coherent mode. Note also that, at relatively low energy transfers, 
H becomes negative, i.e. minimum in June and maximum in December. The negative value, however, for the light 
targets is small for all WIMP masses. For this reason for a light target the integrated modulation amplitude h is 
always positive (maximum in June, minimum in December, as expected). The modulation curves H keep increasing as 
the energy transfer increases, mainly because the time independent amplitude, coming in the denominator, decreases. 
Thus in spite of this increase of H, h remains constant. 



RESULTS FOR THE SPIN CONTRIBUTION 



One amplitude is dominant 



This occurs in cases when the nuclear structure leads to a dominant spin ME, like 19 F with a dominant proton 
component. In this case, barring unusual circumstances at the quark level favoring the component not favored by 
nuclear physics, the analysis is simple. Thus, e.g., in the case of 19 F (f2 p = 1.645, 0„ = —0.030 the event rate for an 
elementary cross section of 10 _5 pb is exhibited as a function of the WIMP mass in Fig. 2] From these plots, for a 
given WIMP mass, one may extract limits on the relevant nucleon cross section from the experimental limits. Using 
the event rate of 13.75 Kg-d or 5020 Kg-y of PICASSO Q and the most favorable WIMP mass of 30 GeV, from 
Fig . [4] we extract a proton spin cross section of 0.1 pb, to be compared with the value of 0.16 pb extracted there 
[24|. From Eq. (f20|) we extract a coherent cross section of 2.5 x 10 _4 pb for this system, which is poor compared to 
the limits of CDMS [6] and XENON Q. The PICASSO people are fighting with their new detector against the a 
background, with a flat plateau in the region of their signal, and their limit will soon substantially improve . 
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5 10 15 20 25 30 

Q -> keV 

FIG. 3: On the left the quantity dRo ^^' m x) an( j on ^he j-jg^ ^hg quantity H(Q,A,m x involving the spin induced 
process for the A=19 system as a function of the energy transfer Q in keV. The thick solid, the dotted , the dashed 
and the thin solid lines correspond to WIMP masses 10,30, 50 and 100 GeV respectively. On the left panel the range 
of Q is restricted to make the curve for low WIMP mass more visible. For masses heavier than 30 GeV the 
differential event rate has essentially a constant slope. So it is adequate to restrict ourselves to low Q. The full 

range of Q can be inferred from the right panel. 




FIG. 4: The event rate (kg-y) for the target 19 F assuming a nucleon cross section of 10 -5 pb as a function of the 
WIMP mass in GeV. In the left panel the continuous curve takes into account only the proton component. The 
dotted curve results when the proton and neutron cross sections are the same, but the corresponding amplitudes are 
opposite (the isoscalar amplitude is assumed to vanish). The difference is small. In the right panel we consider the 
case that the elementary proton cross section vanishes. In this case the nuclear structure suppresses the rate. 



Exclusion plots in the ^/op, ^/on plane 
From the experimental data, using the nuclear spin matrix elements, one can extract a restricted region in the a p , a n 



plane [43j-[44|. The relevant relation is: 



IIOpI^p + |fiJ^ i(W ^| 2 = <ri ~ 3R (25) 

K (<t \)sA 

where sA is a short hand notation for c sp i n (A, /i r (A), m x o). The extracted values, given the event rate and the spin 
ME, depend on the WIMP mass and the relative phase of the two amplitudes. 

Since the procedure is much more complicated than that entering the analysis of the coherent node, a few explanations 
regarding the presentation of our results (Figs [5] fT0|) are in order: 

• We found it more convenient to present in the plots the extracted ^fa^ and ^/a^ rather than the cross sections 
themselves. 

• For illustrative purposes the dependence on 5 can be given in a simple graph whereby the cross sections can be 
expressed in units containing all the parameters. The extracted shapes, which depend on 5a are shown in Fig. 
[5] Such a plot, in principle, contains all the needed information, but it is too general to be practical. 
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(a) (b) 

FIG. 5: A "universal" exclusion plot in the {^/Op, y/on) plane exhibiting the dependence on the phase 8. On tc 
right panel the nuclear spins are of the same sign, while on the other of opposite signs. When the two amplitudes 
are relatively real, they are not bounded except when the the phase 8 is the same with the relative phase of the two 
nuclear matrix elements. In both panels the dotted, the fine solid, the dashed, the dotted- dashed, and the thick 
solid curve correspond to 8 = 0, 7r/6, tt/2, 5 n/6 and n respectively. 



• The contour for 8 ^ 0, tt is in general an ellipse. For a given experimental bound, the allowed values of the 
cross sections are in the space enclosed by an ellipse. One can see that, depending on 8 the maximum allowed 
cross sections can be quite a bit higher than those extracted assuming a single mode. If the two amplitudes are 
relatively real, then the contours become straight lines and the cross sections may be constrained, but only if 
the relative phase of the two amplitudes is the same with 8 a- If they differ by it, the individual cross sections 
are not bounded, they can be anywhere between the two lines. 

• For given nuclear spin ME the extracted Jo~p~ and y/dn are presented in units . / o\ 3J ^ (see Fig. [qI-TTUI ). 

Once the experiment determines the rate R and the parameter sA, for the chosen WIMP mass, is read off from 
table HH one can immediately extract from the figures the cross sections in units of o\ ( K(a\) is given by Eq. 
(|18p). As an illustration we do this on the right panel of the the figures [SI -ITU1 assuming an event rate of 1 event 
per Kg target per year for the optimum value of m x (maximum of sA) 

The following cases are of experimental interest: 

i) We first consider the case of nuclear spin matrix elements of opposite sign and |fi n | > |O p | as is the case of the A=3 
system. The exclusion plots are shown in Fig. [5] 

ii) Next comes the case of spin matrix elements of opposite sign and |f2 p | 3> |f2 n | as is the case of the 19 F target. This 
case has already been analyzed above, considering only protons. Just in case the elementary proton cross section is 
very suppressed, we present the relevant exclusion plots in Fig. [7] 

hi) After this we consider nuclear spin matrix elements of same sign and |O n | 3> |f2 p |. This is the case of the 73 Ge 
target. The relevant exclusion plots are shown in Fig. [5J 

iv We consider the case with both spin matrix elements being significant. Such may be the case of the 127 I target 
(ftp = 1.127, Q„ = 0.315). The resulting exclusion plots are shown in Fig. [9l 

v) We finally show the exclusion plots exhibited by another proton favoring nucleus, 23 Na, which is present together 
with 127 I in the target Nal (see Fig. [T0|) . We notice that, since one has the same number of nuclei of each component 
in a given mass of the target, 23 Na competes well with the 127 I in the spin induced event rate. 



CONCLUDING REMARKS 



We have analyzed the spin induced WIMP nucleus elastic cross section and related event rates. Both depend 
rather sensitively on the spin structure of the nucleus. Barring unusual circumstances at the elementary level, the 
spin mode has no chance to compete with the coherent WIMP nucleus scattering in the case of heavy targets. It 
could, however, compete with it in the case of light targets. For light targets, and in particular for 3 He and 19 F, 
we believe the nuclear matrix elements are very accurate to allow reliable extraction of the nucleon cross sections 
from the data, if and when they become available. In the cases considered here, with the possible exception of 127 I, 



5 10 15 20 

^/dp^ x/1.3 x 10- 4 pb 

(b) 

FIG. 6: The exclusion plot in the (a/cx^, y/on) plane in the case of the target 3 He for various values of the phase 5 
using the relevant spin ME of table IJ The units depend on the parameters of table HH the experimental rate R as 
well as and, for the chosen scale of ci, on K(a%) (a). The same exclusion plot in the case of a WIMP with a mass 20 
GeV normalized to 1 event per kg target per year in the indicated units, obtained using K{u\) = 1.6 x 10 2 y _1 and 
o\ = 10~ 5 pb (b) . Cross sections for other event rates can be trivially extracted by a simple rescaling of panel (b). 
When the two amplitudes are relatively real (5 = 0, it), they are not bounded except when 6 coincides with the 
phase difference 5a of the neclear matrix elements. For the labelling of the curves see Fig: [5] 






(a) (b) 



FIG. 8: The same as in Fig. [6] in the case of 73 Ge target. Now in panel (b) we exhibit the most sensitive case of a 

WIMP mass of 50 GeV. 
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0.2 0.4 0.6 0.8 1.0 1.2 

^a~ v -> Vl-4 X 10- 4 pb 



(a) (b) 

FIG. 9: The same as in Fig. [6] in the case of 127 I target. Now in panel (b) we exhibit the most sensitive case of a 

WIMP mass of 30 GeV. 




the nuclear structure tends to favour the proton or the neutron component. This allows a simple extraction of the 
corresponding nucleon cross section. This is also true even if both components are present, but the isoscalar amplitude 
at the nucleon level is suppressed. Finally, even if both the proton and the neutron amplitudes are important, we 
have shown that knowledge of the nuclear matrix elements allows one to draw suitable exclusion plots. Unfortunately, 
then, the situation is technically a bit complicated by the fact that one must draw one exclusion plot for each WIMP 
mass. So, for targets with spin different from zero, exclusion plots should be drawn as more experimental data become 
available. 
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